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The perspective of the detectabihty of Galactic dark matter subhaloes on the 
Fermi satelhte is investigated in this work. Under the assumptions that dark matter 
annihilation accounts for the "GeV excess" of the Galactic diffuse 7-rays discovered 
by EGRET and the 7-ray flux is dominated by the contribution from subhaloes of 
dark matter, we calculate the expected number of dark matter subhaloes that Fermi 
may detect. We show that Fermi may detect a few tens to several hundred subhaloes 
in 1-year all sky survey. Since EGRET observation is taken as a normalization, this 
prediction is independent of the particle physics property of dark matter. The uncer- 
tainties of the prediction are discussed in detail. We find that the major uncertainty 
comes from the mass function of subhaloes, i.e., whether the subhaloes are "point 
like" (high-mass rich) or "diffuse like" (low-mass rich). Other uncertainties like the 
background estimation and the observational errors will contribute a factor of 2 ~ 3. 

PACS numbers: 95.35.+d 

I. INTRODUCTION 

The dark matter (DM) problem is one of the most important issues in modern physics 
and cosmology. Unfortunately, more than seventy years after the discovery of DM, the 
nature of the DM particle is still unclear. To identify the DM particles, it is necessary to 
"see" them in particle physics experiments beyond the gravitational measurements. There 
are usually three types of experiments suggested to capture the DM particles. The first one 
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is to produce DM particle pairs by high energyparticle colhsion at colhders such as the 



forthcoming Large Hadron Colhder (LHC, 



nergy pan 



The second method is the so-called 



direct search of DM particle by looking for the signals when DM scattering off the detector 



nuclei 



i fl, y, 0^ 



8|, which is thought to be the most direct way to show the existence and 



understand the property of DM particle. Finally the indirect search of DM anni 



products in cosmic rays (CRs) including 7-rays 



neutrinos 



13, 



18 



|9|,[lo|,[ll|,[l 



anti-particles 



illation 
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15| and 



191], is also an important complementary method for DM searches. 



Among the annihilation products, 7-ray is the most attractive one for detection. Com- 
paring with neutrinos, 7-rays are easier to be recorded; while comparing with charged anti- 
particles, 7-rays will not be deflected by the magnetic field and can trace back to the sites 
where the annihilation takes place. In this work we focus on the 7-rays from DM annihila- 
tion. 

It is known that DM annihilation products are proportional to the density square of DM 
distributions. Therefore the highly concentrated regions are good sites for DM searches. 
Theoretically there are many such sites proposed to search for DM, such as the Galactic 
center 
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30 
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23|), DM substructures [24 
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32l |. and the mini-spikes 
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28|, dwarf galaxies [29|, 



351]. It has been recognized that DM subhaloes 



have several advantages for DM searches when comparing with other sites [25|, |36| . Firstly, 
the subhaloes distribute isotropically in the Galactic halo, and can easily be located in a 
low background environment away from the Galactic plane. Secondly, the subhalo may well 
decouple from baryon matter and give clean DM signals. Thirdly, large number of subhaloes 
makes it possible to do statistical study. Therefore to draw attention to the signals from 
DM subhaloes should be a very important aspect for DM indirect searches. 

From the observational point of view, the Energetic Gamma Ray Experiment Telescope 
(EGRET) on board the satellite Compton Gamma Ray Observatory (CGRO) surveyed the 
high energy 7-ray sky in 30 MeV— 30 GeV. EGRET discovered 271 7-ray sources 371] and 



measured the all-sky diffuse 7-ray emission 38 



Among the 271 sources of the third 



EGRET catalog, only 101 ones are identified as known sources. There are 48 additional 
sources that are suggested to be associated with possible counterparts. A substantial fraction 
with over 120 sources is still unidentified. The DM subhaloes are thought to be possibly 
composed of one population of such unidentified sources. However, the analysis of the 
luminosity and spatial distributions show that no more than ~ 20 sources can be attributed 
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to DM clumps 40, |4l|. While for the Galactic diffuse 7-ray emission, the spectrum shows an 



"excess" at energy ~ GeV compared with that predicted by the conventional CR model 

'3, 



38|. 



An appealing interpretation of the "GeV excess" is DM annihilation [42|, |43|, |4J] . By fitting 
the all sky diffuse 7-ray distributions, de Boer et al. found that the all-sky spectra can be 
well reproduced by 50 — 70 GeV supersymmetric DM annihilation, with an isothermal halo 
distribution and two additional ring-like structures 4^1. However, a "boost factor" of ~ 100 
is needed to match the absolute fluxes. The "boost factor" is proposed to be able to come 
from the enhancement of annihilation rate by DM subhaloes. Bi et al. have constructed a 
realistic model to explain both the diffuse 7-ray spectra and the CR observations, taking into 



account the DM subhaloes from numerical simulations [45|, |46|]. It is found that a strongly 
clumpy DM profile is needed to give such a large "boost factor". In order to further test 
this scenario we have to turn to the experiments with higher sensitivity, such as the newly 
launched Fermi satellite^, which will follow the footsteps of CGRO-EGRET to explore the 
high energy 7-ray sky. The Fermi-LAT (Large Area Telescope) is a high quality instrument 
with larger field-of-view (~ 2.5 sr), wider energy bands (10 MeV— 300 GeV) and higher 
sensitivity compared with EGRET. It will be able to improve the sensitivity of indirect 
search of DM annihilation greatly. 

Considering the advantages of detecting DM annihilation from the Galactic subhaloes, 
it is now one of the most prominent scientific goals to search for DM subhaloes by Fermi. 
The perspective of DM detection from subhaloes on the Fermi satellite has been discussed 
extensively in literature (e.g., [471, |48|, |49|, ISOj). Assuming a specific DM model, Baltz et al. 
claimed that for 5-year exposure Fermi might detect ~ 10 DM satellites (subhaloes) with 
significance > 5a 50||. Kuhlen et al. found that based on the numerical simulation of DM 
distributions the number of subhaloes Fermi might detect was from a few to several dozens 
depending on the subhalo boost factor and DM annihilation parameters 49] . Bertone et al. 
showed that DM annihilation around intermediate mass black holes (mini-spikes) might be a 
population of bright 7-ray sources of Fermi 3^ . However, all of these kinds of studies depend 
sensitively on the particle physics model of DM particles and suffer large uncertainties. 

In the present work we try to investigate the performance of DM subhalo detection on 
Fermi, under the assumptions that the EGRET "GeV excess" comes from DM annihila- 
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tion and the DM induced 7-ray emission is dominated by that from the subhaloes. Since 
the EGRET observation of diffuse 7-rays is adopted as a normahzation, this prediction 
is expected to get rid of the uncertainties from the annihilation cross section and 7-ray 
yield spectrum. The degeneracy between the normalization factor of DM subhaloes and 
the annihilation cross section is discussed. Furthermore we discuss the uncertainties of this 
prediction. 

The outline of the paper is as follows. We first introduce the model of DM annihilation 
scenario to explain the EGRET "GeV excess" in Sec. 2. In Sec. 3 we give the prediction 
of detectable number of subhaloes on Fermi based on this model. The uncertainties of the 
prediction are carefully discussed in Sec. 4. Finally we give conclusion and discussion in 
Sec. 5. 



II. THE MODEL TO ACCOUNT FOR THE EGRET OBSERVATIONS 



One way to explain the EGRET "GeV excess" is the DM annihilation scenario 42l. l43l. l44 1 



This scenario has been extended to be a more realistic one to explain both the diffuse 7-ray 
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461]. For the convenience of discussion in this work, we base on 



and CR observations 

the model of Ref.j46[, however, the conclusion will not be limited in this detailed model. In 
the model of Ref . [46| , the background 7-rays are calculated using the package GALPROP 



51 



52l |. The parameters of GALPROP are adjusted to reproduce the EGRET data and the 



CR observations, especially the antiproton data, after including the contribution from DM 
annihilation. In Fig{T]we show the differences between the observational spectra by EGRET 
and the predicted background (including the extra-galactic component) at different sky 
regions. The six sky regions are defined following Ref. 53|]: region A corresponds to the 
inner Galaxy with 0° < / < 30°, 330° < / < 360° and |6| < 5°; region B is the Galactic 
plane excluding region A with 30° < / < 330° and \b\ < 5°; C is the outer Galaxy with 
90° < / < 270° and \b\ < 10°; regions D, E, F cover the medium and high latitudes at all 
longitudes with 10° < |6| < 20°, 20° < |6| < 60° and 60° < \b\ < 90° respectively. From 
the figure we notice that the differences of the spectra at different regions are very similar, 



which implies that t 



DM annihilation 



lere may be some common origins of the lacked 7-rays, for example the 



42 



43 



44|. 



The diffuse 7-ray fluxes from DM annihilation are determined by two factors: the "particle 
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FIG. 1: The differences between the observational EGRET fluxes and the CR induced 7-rays 
predicted from GALPROP. Similar shapes of different sky regions indicate a common origin of 
the excess. The solid line shows the DM annihilation spectrum (with arbitrary normalization) for 
neutralino with mass 50 GeV. 



factor" which describes the particle physics property of DM particles and the "astrophysics 



factor" which describes the spatial distribution of DM. In Ref. j46| the "particle factor" is 
calculated under the framework of the minimal supersymmetric extension of the standard 
model (MSSM, 5J]). The lightest particle neutralino in MSSM is adopted as the DM particle 



^. The annihilation cross section and 7-ray production spectrum are calculated using Dark- 
SUSY ^ a package for MSSM DM calculation 



A random scan in the MSSM parameter 



space is performed and the adopted models are chosen to satisfy the relic density constraint 



and give large 7-ray fluxes. The annihilation spectrum for the neutralino with mass 



50 



GeV, as shown in FigH], is found to be in good agreement with the observations. 

The "astrophysics factor" is mainly determined according to the numerical simulation of 
DM distribution. High resolution simulations show that a fraction of 10% ~ 20% of the total 



^ However, it will be shown in the following that this assumption can be relaxed in the present work. 
^ http://www.physto.se/ edsjo/darksusy/ 
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mass survives in the self-bound substructures [56|, [STL |58|, l59|, l60|, iMi- The substructures of 



DM can effectively enhance the annihilation signals [2j, |27| . The number density distribution 
of subhaloes from simulation can be fitted by an isothermal spatial distribution and a power- 
law mass function 

iv„f??^V°— (1) 



dN 



dmsub ■ 4:TTr^dr 



\ Mvir / 1 + (r/rn) 



where Mvir ~ 10 M0 is the virial mass of the Galaxy, rn ~ 0.14rvir = 29kpc is the core 



radius of the distribution of subhaloes 



62l ]. and Nq is the normalization factor determined 



by the total mass of subhaloes. The slope a varies from 1.7 to 2.1 in various works 57|, |58|, 



60 



63 



64 



65l |. The intermediate value a ~ 1.9 is favored by the recent highest resolution 



simulation Via Lactea 6^. In Ref. 46|] a = 1.9 is adopted. The masses of subhaloes range 
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67 



to 



from the minimum ~ 10~ M0 which is close to the free-streaming mass 
the maximum about O.OlMyir. The DM density distribution inside the subhalo follows a 7 
profile p(r) = ^^f^f^^.,^. y^^^ with 7 = 1.7, which is steeper in the center than the NFW 



type {(X r ^, [69|]) and Moore type (oc r 



-1.5 



TOQ. A steeper inner slope for small mass 



halo is suggested by numerical simulations 7l|, l72]. The profile parameters rg and Ps are 
determined according to the halo mass M^ir and the concentration parameter Cvir jisl, S] • 



731 ]. Finally the central 



The Cvir — M^ir relation is adopted as the model of Bullock et al. 
density of a subhalo is truncated at a certain radius to avoid the divergence. A characteristic 
radius rcut is introduced within which the DM density is kept a constant Pmax due to the 
balance between the annihilation rate and the in-falling rate of DM Typically we have 

Pmax = lO^S ~ 10^9 kpC-3 H- 

This configuration of DM distribution is then used to calculate the diffuse 7-ray emission. 
Taking Pmax = 10^^ Mq kpc"^ and subhalo mass fraction 20% ^ we find that the results can 
well reproduce the high latitude observations by EGRET (i.e., F region). It can be noted 
from Fig. 2 of Ref. j46| that for the direction far away from the Galactic center the contribution 
to 7-rays from subhaloes dominates that from the smooth halo. For example at the Galactic 
pole direction b = ±90°, the ratio of the astrophysics factors between the subhaloes and 
smooth halo is / = $sub/'^'sm ~ 10^, which is consistent with the requirement of a boost 



factor ~ 100 of Ref. 



42j. For other sky regions especially for the Galactic plane re gion s, the 



enhancement of subhaloes is not enough, and two additional DM rings are needed 



4^. The 



Note that the parameter adoption is a bit different from Ref. [46| , where r^ut instead of Pmax is used 
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overall results including the DM rings are well consistent with the observations. 

III. DETECTION PERFORMANCE OF SUBHALOES ON FERMI 

In this section we calculate the 7-ray fluxes from massive DM subhaloes adopting the 
same DM model presented in Sec. 2 and discuss the detectability of these 7-ray sources on 
Fermi. The Monte— Carlo (MC) realization method is adopted to generate the DM subhaloes 
with mass > 10^ Mq in the Milky Way, following the distribution function Eq.([l]). There 
are about 1.5 x 10^ subhaloes found with mass heavier than 10^ Mq. Totally 100 Milky 
Way like galaxies are generated. For each realization, we calculate the annihilation flux of 
each subhalo with energy threshold -Eth = 100 MeV and count the accumulative number as 
a function of the threshold flux $. The result is shown in Figl2j 

The sensitivities of EGRET and Fermi for 1-year all-sky survey at 5a are 5 x 10~® ph 
cm-2 s"^ (> 100 MeV) Q and 4 x 10"^ ph cm'^ s"^ (> 100 MeV) respectively^ which are 
shown by the vertical lines in FiglS From the flgure we flnd that the number of subhaloes 

4o|,|41|. 



that EGRET can detect is 18.7 ± 4.4, which is consistent with the results of Refs. 
For Fermi, the detectable number is 245.2 ± 16.8, which is an order of magnitude more than 
EGRET. The scattering comes from different realizations. The probability distribution 
of the detectable numbers on Fermi is well fltted with a Gaussian distribution, as shown 
in the top-left panel of FigOl We also plot the mass (log^ol^^sub/^o]) distribution, the 
distance distribution and the directional skymap in Galactic coordinate of the detectable 
DM subhaloes on Fermi for one of the realizations in FigjSl It is shown that subhaloes 
with masses ~ 10® Mq and distances within ~ 50 kpc are more likely to be detected by 
Fermi. The direction distribution is isotropic, which will be signiflcantly different from other 
astrophysical populations of sources. 

In Figl2]the threshold energy of the detector is adopted as 100 MeV. However, since the 
energy spectrum of DM annihilation is known, as shown in Fig{T], the detector performance 
can be optimized by taking a proper energy cut. In FigJH we show the ratio between the 
integral flux of DM annihilation spectrum and the integral sensitivity of Fermi. It will be 
more efficient for the DM detection for larger ratio. It can be seen that for the energy 
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FIG. 2: The accumulative number of subhaloes as a function of integral flux for energy threshold 
E'th = 100 MeV. The two vertical lines show the sensitivities of EGRET and Fermi. 



threshold Eth ~ 1 GeV, the detectability is most optimized. We show that for Eth = 1 GeV 
the number of subhaloes that can be detected by Fermi is 620.4 ± 23.1. In the following of 
this paper we will adopt the optimized threshold energy Et\^ = 1 GeV for discussion. 



IV. UNCERTAINTIES 



In the previous section, we show our prediction of the detectable number of DM subhaloes 
on Fermi based on the EGRET observations of the Galactic diffuse 7-r ays. It is about 250 
(or 620 for optimized energy cut) according to the model of Ref. 46|. In this section we 
discuss the possible uncertainties of this prediction. 



A. Particle factor 



There are two quantities in the particle factor to affect the annihilation flux: the velocity 
weighted averaged cross section {av) and the 7-ray production spectrum per annihilation 
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FIG. 3: Distributions of the detectable DM subhaloes on Fermi: detectable number distribution due 
to various realizations (top-left); mass (top-right), distance (bottom-left) distributions and skymap 
(bottom-right) of the detectable DM subhaloes in one realization. The area of the circle in the 
skymap is proportional to the flux. 



dN/dE. In Refs. 42, |46|], the neutralino is adopted as the DM particle, (av) and dN/dE 



are calculated under the MSSM model. However, from Fig{T]it can be seen that the annihi- 
lation spectrum is fixed by the "gaps" between the observations and the background, and is 
independent with the DM particle models. We find that if we adopt a model independent 
spectrum by averaging the spectra of the 6 sky regions (normalized by the maximum value 
and extrapolate logarithmically to higher and lower energies) instead of the spectrum from 
neutralino annihilation, the results are almost the same. 

The annihilation cross section (^av) will affect the absolute flux of 7-rays. Since the 
absolute flux is determined by the EGRET observation, there is degeneracy between the 
cross section and the astrophysics factor. That is to say, a smaller cross section (erf) needs 
to be compensated by a larger astrophysics factor. Therefore, if the total number of DM 
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FIG. 4: Integral signal from DM annihilation to sensitivity of Fermi ratio as a function of threshold 
energy. 

subhaloes maintains the same, varying (av) is not expected to change the detectabihty of 
DM subhaloes significantly. We will show in the next section that even the normalization 
of the total number changes, the conclusion still holds. However, it should be pointed out 
that, if (av) is so large that the smooth contribution dominates the 7-ray emission from DM 
annihilation (i.e., / = $sub/'^'sm ^ 1); it "will be more difficult to detect DM subhaloes. 

B. Astrophysics factor 

In the astrophysics factor, the normalization of number of subhaloes, the mass function of 
subhalo number distribution, the inner profile of DM subhalo and the concentration model 
may all lead to uncertainties of the predicted detectable number of subhaloes given in Sec. 
III. We now investigate these issues one by one in detail. 

The normalization of the total number of subhaloes, i.e., Nq of Eq.([T]), is determined by 
the mass fraction of clumps in the Galactic halo. It shows relatively large uncertainty in 
different simulations. However, we find that even though the factor A'o is left free in the 
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range of one or two orders of magnitude, the result of this work is still unchanged. The 
reason is that in order to keep the 7-ray emission of region F unchanged, varying A^^q by 
some factor needs to be compensated by the same factor of the annihilation flux of each 
subhalo (e.g., through the rescale of (erf)). From Figj2] we can see that the cumulative 
number A^(> $) oc which means that the change of A^o is equivalent to shifting the 
curves in Figl2] downward (or upward) while the change of (av) shifts the curves rightward 
(or leftward) by the same factor. Therefore the number of subhaloes with flux higher than 
the sensitivity of Fermi will keep the same. 

The index of the mass function a in Eq. ([T]) will affect the mass distribution of subhaloes, 
and accordingly affect the A^(> $) ~ $ relation. If a is smaller, the fraction of contribu- 
tion to the diffuse 7-rays from high mass subhaloes becomes more important (see Fig. 4 of 
Ref . {tgI ) . so we can expect that the subhaloes with higher fluxes are richer. The expected 
numbers of subhaloes on Fermi for energy threshold 1 GeV, for different a are listed in Table 
m For each case, we scale (av) to recover the EGRET observation at region F. The inner 
DM profile 7 = 1.7, the central maximum density pmax = iC^^ Mq kpc~^ and subhalo mass 
fraction 20% are kept unchanged ^. It is shown that for different values of a the results 
differ from each other significantly. The parameter a determines the DM to be high mass 
rich ("point like") or low mass rich ("diffuse like"), which is the key factor for the detection 
of DM subhaloes. It should be noted that for the cases a = 1.7 and 1.8, the most luminous 



DM subhalo is even brighter than the brightest unidentified EGRET source STj. In order 
not to break the constraint from the EGRET source catalog we find a > 1.9. The index of 
the mass function of subhaloes a is actually the most important source of the uncertainties 
when predicting the detectability of DM annihilation from subhaloes. 

TABLE I: Detectable number of DM subhaloes on Fermi for energy threshold 1 GeV. 



a 1.7 


1.8 


1.9 


2.0 


2.1 


Number 4068 ± 32.8 


2147 ± 40 


620.4 ±23.1 


110.9 ± 10.5 


13.6 ±3.6 



The inner slope of DM subhalo can also result in slightly different mass dependence of the 



^ The same treatment is also employed in the following discussion, i.e., we rescale (av) to cancel the effect 
by the relevant changes and keep other settings unchanged as the benchmark model presented in Sec. 2. 
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annihilation flux, which may also change the detection performance of subhaloes on Fermi. 
This effect is not very important. We show that for Moore profile, iV^et = 410.8±20.0, and for 
NFW profile, iVdet = 441.4 ± 20.3. From Fig.3 of Ref.0 we can see that for 7 = 1.7 profile 
the fraction of DM annihilation luminosity from high mass subhaloes is more important than 
that of Moore or NFW profiles. Therefore the detectable number 620.4 ± 23.1 for 7 = 1.7 
profile (i.e., the benchmark model in this work), is some larger than that of NFW or Moore. 



Finally, the concentration model adopted in Ref . |46l] is from Ref . [731] . At present there is 
no good determination about the subhalo concentration, especially for the low mass haloes 
which are beyond the resolution of simulations. As a comparison, using another concentra- 
tion model of Eke et al. 781] we find A'dot = 794.6 ± 26.3. The model of Eke et al. 78|] gives 
larger concentration (accordingly larger annihilation flux) of heavier halo, resulting in more 
detectable number of subhaloes. 



C. Background estimation 

The background 7-ray emission from CRs is calculated based on the conventional CR 
propagation model. To discuss the uncertainty of the background calculation is very com- 
plicated due to the uncertainties of the inputs, including gas distribution in the Galaxy, 
interstellar radiation field, nuclear interaction cross section, CR measurements, CR source 



distribution, solar modulation and so on 53|]. However, many of these factors are de 



generated. For example a lower gas density can be compensated by a higher CR source 
normalization. According to the CR and diffuse 7-ray measurements, the uncertainty of 
the estimation of 7-ray background can be limited to a few tens percent. Here we adopt a 
rough estimation of about 20 ~ 30% uncertainty on background 7-rays based on the model 



fitting of propagation parameters using B/C ratio data 79|]. If we change the background 
intensity by ±25%, the corresponding DM contribution will change ~ ±50% to compensate 
this variation. The detectable number of DM subhaloes will also change by ±50%. 

D. Observational uncertainties 



The systematic errors of EGRET are estimated as ~ 15%. If the statistical errors of 
the photon counts are taken into account the errors are larger. It is shown that even 
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the DM contribution varies ±50%, the total 7-ray emission is still consistent with EGRET 
observation in la level. This will correspond to ~ 50% uncertainty of the detectable number 
of DM subhaloes. 

V. CONCLUSION AND DISCUSSION 

In this work we study the detection performance of DM subhaloes on Fermi under the 
assumption that the "GeV excess" discovered by EGRET comes from the DM annihilation. 
Considering that EGRET being the precursor of Fermi it is quite natural to predict the 
perspective of Fermi based on the EGRET results. An additional assumption is that the 
exceeded part of diffuse 7-rays at high Galactic latitude is dominated by the contribution 
from DM subhaloes, i.e., DM subhaloes give a large boost factor. Then the EGRET ob- 
servation is used as a normalization and we get the detection perspective of DM subhaloes 
on Fermi. We find the possible detectable number of DM subhaloes is from a few tens to 
several hundred after 1-year operation of Fermi. For the favored value of the mass function 
slope of DM subhaloes a ~ 1.9 from recent highest resolution numerical simulation, there 
will be several hundred subhaloes that can be seen. 

This prediction has relatively large uncertainties. A major uncertainty comes from the 
mass function of DM subhaloes, which is not well determined by numerical simulations. 
The mass function determines the ratio of 7-ray fluxes between low mass subhaloes and 
high mass ones. If the mass function is steep, i.e., low mass subhaloes dominate the 7-ray 
emission, then the 7-ray sky will be more "diffuse like" and the detectable number of DM 
point sources decreases. Otherwise more DM subhaloes can be detected. Our result shows 
that according to the current knowledge about the mass function from numerical simulations 
the detectable number can vary by orders of magnitude. Other uncertainties from the inner 
properties of DM subhalo, background estimation and the observational errors will totally 
contribute a factor of 2 ~ 3. Different from previous studies on the DM indirect search, our 
result is independent of the particle physics property of DM particles. 

In spite of the large uncertainties, our result shows that searching for DM subhaloes may 
be a promising way for the indirect search of DM on Fermi. It is believed that Fermi will 
open a new era in the DM study and greatly enrich our knowledge about DM. 

Finally we point out that this prediction is based on the assumption that the EGRET 
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le CR processes contribute 



53| , or the "GeV excess" is 



"GeV excess" completely comes from the DM annihilation. If t 
a fraction or the whole of the "GeV excess" as proposed in Ref. 
due to the wrong sensitivity estimation of the detector of EGRET as pointed out in Ref. 80| , 
the detectability of DM subhaloes on Fermi will not be so promising as we have discussed. 
Especially the prediction becomes difficult due to the large uncertainties from the particle 
physics and distribution of DM. 
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